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ABSTRACT 

We investigate gas metallicity of narrow-line regions in narrow-line Seyfert 1 galaxies (NLSls) and 
broad-line ones (BLSls) in order to examine whether or not there is a difference in the gas metallicity 
between the two populations of Seyfert 1 galaxies. We apply two methods to study this issue. One 
is to use the emission-line flux ratio of [N n] A6583/Ha narrow in combination with some other optical 
emission-line flux ratios. This method, which has been often applied to Seyfert 2 galaxies, suggests 
that the gas metallicity of narrow-line regions is indistinguishable or possibly higher in BLSls than in 
NLSls. On the contrary, the other method in which only forbidden emission-line fluxes are used results 
in that NLSls tend to possess metal-richer gas in the narrow-line regions than BLSls. We point out 
that this inconsistency may be owing to the contamination of the broad component of permitted lines 
into the narrow component of ones in the first method. Since the results derived by using only forbidden 
emission-line fluxes do not suffer from any uncertainty of the fitting function for the broad component 
of Balmer lines, the results from this method are more reliable than those derived by using permitted 
lines. We thus conclude that the gas metallicity of narrow-line regions tends to be higher in NLSls than 
in BLSls. 

Subject headings: galaxies: abundances - galaxies: active - galaxies: nuclei - galaxies: Seyfert - 
quasars: emission lines 



1. INTRODUCTION 

Seyfert nuclei are typical active galactic nuclei (AGNs) 
in the nearby Universe. They have been broadly classified 
into two types based on the presence or absence of broad 
(;> 2000 km -1 ) permitted lines in their optical spectra 
(Khachikian & Weedman 1974); Seyfert nuclei with broad 
lines are type 1 (hereafter SI) while those without broad 
lines are type 2 (S2). This difference can be explained 
by the "unified model" in which it is considered that the 
broad-line region (BLR) is surrounded by a dusty torus 
and that the BLR in S2s is obscured by the edge-on view 
torus from the observer (e.g., Antonucci & Miller 1985; 
Antonucci 1993). Therefore only narrow (<; 2000 km -1 ) 
lines, which arise from narrow- line regions (NLRs), are 
seen in the spectra of S2s. In addition to these two types of 
Seyfert nuclei, narrow-line Seyfert 1 galaxies (NLSls) also 
comprise a sub-class of the Seyfert nuclei. Since NLSls 
exhibit only narrow (i.e., < 2000 km -1 ) permitted lines in 
their optical spectra, NLSls are recognized as a distinct 
type of ordinary broad-line Sis (BLSls). However, NLSls 
also show several characteristics of Sis (e.g., Osterbrock & 
Pogge 1985); e.g., small line flux ratios of [O m]A5007/H/3 
(<; 3), strong Fe n and some high- ionization emission lines, 
and non-absorbed X-ray spectra. Therefore, it is believed 
that the BLR of NLSls is not obscured but the line width 
of the BLR is narrow for some reason. 

One possible idea to explain the origin of the narrow 
BLR emission of NLSls is that NLSls have a smaller mass 
black hole than BLSls (e.g., Boiler, Brandt, & Fink 1996; 
Laor et al. 1997). If the BLR in NLSls is the same as 
the BLR in BLSls but the gas motion is more quiescent 
because of the smaller black hole mass, the narrowness 



of the BLR emission can be successfully explained. The 
steep and highly variable X-ray spectra, which are also 
the common properties of NLSls (e.g., Boiler et al. 1996; 
Turner et al. 1999a; Leighly 1999a, 1999b), can be also 
explained by introducing the small black hole mass (e.g., 
Wang, Brinkmann, & Bergeron 1996; Hayashida 2000; Mi- 
neshige et al. 2000; Lu & Yu 2001; Puchnarewicz et al. 
2001). 

Recently, it has been proposed that NLSls may be rela- 
tively young AGNs with a black hole still in a growing-up 
phase (e.g., Mathur 2000a, 2000b; Mathur, Kuraszkiewicz, 
& Czerny 2001; see also Komossa & Mathur 2001; Mathur 
2001). If NLSls possess a smaller mass black hole than 
BLSls while the luminosity is similar between NLSls and 
BLSls, a mass accretion rate for the same black hole mass 
is larger in NLSls than in BLSls. This implies that the 
mass accretion in NLSls is more efficient than in BLSls, 
which leads to the above hypothesis. Indeed some observa- 
tions suggest that the mass accretion rate in NLSls is high, 
close to the Eddington rate (e.g., Pounds, Done, & Os- 
borne 1995; Kuraszkiewicz et al. 2000; Pounds & Vaughan 
2000; Wandel 2000; Puchnarewicz et al. 2001). Moreover, 
Wandel (2002) reported that NLSls have the smaller black 
hole mass than the value inferred from the correlation be- 
tween the black hole mass and spheroidal mass of the host 
galaxy (see also Wandel 1999, 2000; Mathur et al. 2001), 
which is consistent with the hypothesis that NLSls possess 
a growing black hole in their nucleus. 

In the framework of this scheme, the gas metallicity 
of NLSls has been discussed very recently. Wills et al. 
(1999) found a negative correlation between the emission- 
line width of H/3 and the line flux ratio of N v A1240/C 
iii]A1909, which can be interpreted as a result of a nitrogen 
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Table 1 

Statistical Properties of Narrow Emission-Line Flux Ratios' 1 of the NLSls and the BLSls 



Line Ratio 


NLSl b 

average 


median 


average 


BLSl b 

median 


Pks 


[N n]A6583/Ha 

[0 i]A6300/[N n]A6583 

[0 n]A3727/[N n]A6583 

[S n]AA6717,6731/[N n]A6583 


0.406 
0.123 
0.281 
0.427 


0.234 
0.100 
0.237 
0.423 


0.462 
0.261 
0.686 
0.590 


0.367 
0.204 
0.489 
0.433 


0.912 
0.133 
0.533 
0.767 


a Dust extinction is not corrected 












b Data are taken from Rodrfguez- 


Ardila et al. (2000) 











Table 2 

Statistical Properties of the Estimated Oxygen Abundances of NLRs in the NLSls and BLSls 





average 


NLS1 


median 


average 


BLS1 


median 


12 + log (0/H)sscKi a 
12 + log (0/H) SSCK2 a 


8.420 
8.679 




8.384 
8.739 


8.517 
8.654 




8.512 
8.752 



a Calculated by equations (1) and (2), which are derived by Storchi-Bergmann et al. (1998). 



overabundance in BLRs of NLSls (see also Wills, Shang, 
& Yuan 2000; Shemmer & Netzer 2002). The observed 
strong fluorescent Fe Ka and optical Fe n emission, and 
the unusual spectral feature at 1 keV, which are seen in 
spectra of some NLSls, may suggest the overabundance of 
iron in NLSls (Turner, George, & Netzer 1999; Ulrich et 
al. 1999; Collin & Joly 2000; see also Boiler et al. 2002). 
In order to investigate the gas metallicity of NLSls and 
BLSls further, it seems crucial to examine the gas metal- 
licity in NLSls via various approaches. Therefore, in this 
paper, we study the gas metallicity of NLRs in NLSls and 
compare it with that in BLSls using some optical diagnos- 
tic emission-line flux ratios. 

2. DATA 



reason for this is that the accurate measurement of the 
narrow component of permitted lines is rather difficult for 
Sis. Thus high-quality spectra and their careful analy- 
sis are required to measure the narrow-component flux of 
permitted lines. 

Recently, Rodrfguez-Ardila, Pastoriza, & Donzclli 
(2000) presented high-quality optical spectra of 7 NLSls 1 
and 16 BLSls 2 . They carefully removed the underlying 
stellar continuum and multiple Fe n emission lines from 
the observed spectra. Then they deconvolved the narrow 
component of permitted lines from the broad component 
and measured various emission-line fluxes accurately. We 
basically refer the data 3 of Rodriguez- Ardila et al. (2000) 
in this paper to investigate the gas metallicity of NLRs in 
NLSls and BLSls. 

In order to see some possible correlations among narrow 
emission-line flux ratios in a larger sample of Seyfert nu- 
clei, we also refer the data compiled by Nagao, Murayama, 
& Taniguchi (2001c). Their sample contains «350 Seyfert 
nuclei. Since their work is based on the data compilation 
from the literature, the measurement of emission-line flux 



The gas metallicity in NLRs of AGNs has been often 
investigated for S2s (e.g., Shields & Oke 1975; Storchi- 
Bergmann & Pastoriza 1989, 1990; Storchi-Bergmann 
1991; Alloin et al. 1992; Kraemer et al. 1994; Schmitt, 
Storchi-Bergmann, & Baldwin 1994; Radovich & Rafanclli 
1996; Kraemer, Ruiz, & Crenshaw 1998) while it has been 
scarcely studied for Sis (e.g., Kraemer et al. 1999). One 

!CTS H34.06, Mrk 1239, CTS J03.19, CTS J04.08, NGC 4748, CTS J13.12, and 1H 1934-063. 

2 CTS C16.16, MCG -5-13-17, CTS H34.03, CTS B31.01, Fairall 1146, CTS M02.30, CTS J07.02, CTS J10.09, CTS R12.15, CTS J14.05, 
CTS J15.22, CTS M17.17, CTS G03.04, 1H 2107-097, CTS A08.12, and CTS F10.01. Note that we do not use the data of CTS A08.12 in the 
following analysis since the flux of the [N n]A6583 emission of this object is not presented in Rodrfguez-Ardila et al. (2000). 

3 Note that there is a type-setting error in the Rodrfguez-Ardila et al. (2000) paper; i.e., a second header of Table 4 in the Rodrfguez-Ardila 
et al. (2000) paper was cut in the editing process. The second set of line fluxes (when the object names are repeated) corresponds in fact to 
other emission lines organized as follows: [O i]A6300, [Fe x]A6374, H aA6563, [N n]A6584, [S n]A6717, [S n]A6731, [Ar in]A7135, [O n]A7325, [S 
in]A9069, and E(B - V). 
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Fig. 1. — Frequency distributions of the emission-line flux ratio of [N ll]A6583/Ha for the NLSls and the BLSls in the sample of 
Rodriguez- Ardila et al. (2000). The data are not corrected for the dust extinction. 



ratios is done in heterogeneous ways. Therefore, we refer 
their results only when we see general trends of emission- 
line flux ratios and we never use them to discuss the differ- 
ence in the properties of the NLR emission between NLSls 
and BLSls. 

3. PREVIOUS DIAGNOSTICS FOR THE GAS METALLICITY 

OF NLR 

3.1. The [Nn]\6583/Ha Method 

The gas metallicity in NLRs of S2s has been some- 
times studied by using the emission-line flux ratio of [N 
n]A6583/Ha (e.g., Storchi-Bcrgmann & Pastoriza 1989, 
1990; Storchi-Bcrgmann 1991; Radovich & Rafanclli 
1996). Using photoionization models, Storchi-Bcrgmann 
& Pastoriza (1989, 1990) showed that this emission-line 
flux ratio is sensitive to the nitrogen abundance but rather 
insensitive to the gas density and the ionization parame- 
ter (i.e., a number density ratio of ionizing photons to 
hydrogen atoms). Note that the nitrogen abundance is 
important since the N/O abundance ratio in galactic H II 
regions is known to scale with the O/H abundance ratio 
(e.g., Shields 1976; Pagel & Edmunds 1981; Villa-Costas 
& Edmunds 1993; van Zee, Salzer & Haynes 1998; Izo- 
tov & Thuan 1999). Therefore the nitrogen abundance, 
N/H, scales with Z 2 , where Z is the gas metallicity. The 
data of H n regions show that these scaling relations hold 
approximately for Z J> 0.2Z Q . 

We show the histograms of the emission-line flux ra- 
tio of [N ii]A6583/Ha narrow for the NLSls and the BLSls 
in the sample of Rodriguez- Ardila et al. (2000) in Fig- 
ure 1. The data arc not corrected for the dust extinction 
because the effect of the dust extinction on this flux ra- 
tio is negligibly small. The average and median values 
of [N ii]A6583/Ha narrow for the NLSls and those for the 
BLSls are given in Table 1. There is no apparent differ- 
ence in this flux ratio between the two populations. It 
is also suggested by the Kolmogorov-Smirnov (KS) sta- 
tistical test that the two frequency distributions of the 
flux ratio of [N n] A6583/Ha narrow are statistically indis- 



tinguishable (Pks = 0.912). As shown in Figure 2, we 
cannot find an apparent correlation between the flux ra- 
tio of [N ll]A6583/Ha narrow and full- width at half maxi- 
mum (FWHM) of the broad component of the Ha emis- 
sion. The corresponding correlation coefficient is 0.009, 
which implies that there is no meaningful correlation be- 
tween the two quantities. This result seems contrary to 
the result of Wills et al. (1999), who found an apparent 
negative correlation between the emission-line flux ratio of 
N v A1240/C in] A1909, which is sensitive to the nitrogen 
abundance, and the FWHM of the broad component of 
the H/3 emission. 

3.2. The SSCK Method 

Based on the intensive photoionization model calcu- 
lations and the careful calibrations with observations, 
Storchi-Bergmann et al. (1998, hereafter SSCK) proposed 
the revised diagnostics for the gas metallicity of NLRs in 
AGNs. They showed that the oxygen abundance of the 
NLR gas in AGNs can be expressed by 

12 + log(0/H) SSC Ki = 8.34 + 0.212a; - 0.012a; 2 

- 0.002y + 0.007a;?/ - 0.002a; 2 y 

+ 6.52 x 10"V 

+ 2.27 x 10~ 4 a;y 2 

+ 8.87 x lO^'Vy 2 (1) 

where x = [N n]A6583/Ha narrow and y = [O 
lll]A5007/H/3 narrow , and 

12 + log(0/H)gscK2 = 8.643 - 0.275u + 0.164u 2 
+ 0.655t;- 0.154to 
- 0.02lu 2 v + 0.288i; 2 
+ 0.162m; 2 + 0.0353uV (2) 

where u = log ([O n]A3727/[0 m]A5007) and v = log ([N 
ii]A6583/Ha na rrow)- These methods weakly depend on the 
gas density; the derived values should be subtracted by 
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Fig. 2. — Diagram of the flux ratio of [N n]A6583/H« versus the FWHM of the broad component of the H« emission. The circles and the 
triangles denote the NLSls and the BLSls of the sample of Rodrigucz-Ardila ct al. (2000), respectively. 



the correction term, 0.1 x log (nu/300 cm~ 3 ) where njj 
is the hydrogen density (see SSCK for more details). We 
apply these relations to the data of Rodriguez- Ardila et al. 

(2000) . For the objects whose gas density is not given in 
Rodriguez- Ardila et al. (2000) (CTS H34.03, CTS M02.30, 
CTS J14.05, CTS G03.04, and 1H 2107-097), we assume 
the gas density to be 300 cm~ 3 . The average and me- 
dian values of the derived 12 + log (0/H)sscki and 12 + 
log (0/H)sscK2 are given in Table 2. Here the effect of 
the dust extinction is corrected by adopting the values of 
Ay presented by Rodrigucz-Ardila ct al. (2000) and the 
extinction curve of Cardelli, Clayton, & Mathis (1989). 
There is no tendency that the NLSls are metal rich com- 
pared to the BLS1 as implied by some earlier studies. To 
the contrary, the derived oxygen abundances may be larger 
in the sample of BLSls than in the sample of NLSls. We 
show the relationships between the FWHM of the broad 
component of the Ha emission and the derived oxygen 
abundances in Figure 3. The object with a larger FWHM 
of the Ha emission seems to be more oxygen abundant 
when equation (1) is applied, though this tendency is not 
significant. On the other hand, such a tendency is not seen 
when equation (2) is applied. 

3.3. Remarks on the Two Methods 

Do the above results presented in §3.2 suggest that 
the NLR of NLSls is not metal rich compared to that 
of BLSls? Recently, Veron-Cetty, Vcron, & Gongalves 

(2001) pointed out that the broad component of permit- 
ted lines in the optical spectra of NLSls is well fitted by a 
Lorentzian profile and thus the decomposition of the nar- 
row component of permitted lines from the broad compo- 
nent by Rodrigucz-Ardila et al. (2000) may be inaccurate, 
since the broad components of permitted lines were fitted 
by a Gaussian profile by Rodriguez-Ardila et al. (2000) 
(see also Moran, Halpern, & Hclfand 1996; Leighly et al. 



1999b; Goncalves, Veron, & Veron-Cetty 1999; Sulentic 
et al. 2002). If this claim is the case, the flux ratio of [N 
ii]A6583/Ha narrow may be inappropriate for the metallicity 
diagnostics because the selection of the fitting function for 
the broad component of permitted lines affects the mea- 
surement of not only the broad component but also the 
narrow component of permitted lines. 

Indeed, the derived oxygen abundances of the NLR gas 
of the NLSls and the BLSls by using equations (1) and (2) 
should be suspected if we take account of the following two 
facts. First, most of the derived oxygen abundances sug- 
gest very sub-solar metallicities. This is especially remark- 
able when equation (1) is adopted: the median value of the 
derived oxygen abundances corresponds to Z ~ 0.32Z & . 
Second, the derived oxygen abundances are not consistent 
to each other. That is, the average and median values of 
the derived oxygen abundances are systematically differ- 
ent: the oxygen abundances derived by using equation (2) 
are larger than that by using equation (1) by ~0.3 dex. 
This difference is more clearly shown in Figure 4, in which 
we compare the two kinds of the inferred oxygen abun- 
dances. Although SSCK also mentioned that the oxygen 
abundance derived by equation (2) tends to be higher than 
the one derived by equation (1) by ~0.11 dex, the differ- 
ence in the sample of Rodriguez-Ardila et al. (2000) is 
much larger than that reported by SSCK. This may imply 
that the methods proposed by SSCK are inappropriate to 
estimate the gas metallicities of NLRs, at least for Sis. 
Thus it seems to be safe to avoid using the narrow com- 
ponent of permitted lines for the investigation of the gas 
properties of NLRs in Sis. 

4. NEW DIAGNOSTICS FOR THE GAS METALLICITY OF 

NLR 

4.1. Correlations among Forbidden Emission-Line Flux 

Ratios 
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Fig. 3. — Diagram of the oxygen abundances estimated by the equations given by Storchi-Bcrgmann ct al. (1998) versus the FWHM of 
the broad component of the Ha emission. The symbols are the same as those in Figure 2. 
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Fig. 4. — Diagram of the oxygen abundance estimated by the equation (1) versus that by the equation (2). The loci of equal values for 
the two quantities are drawn by a dot line for comparison. The symbols are the same as those in Figure 2. 
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Fig. 5.— Diagrams of emission-line flux ratios of (a) [O n]A3727/[N n]A6583 versus [O i]A6300/[N n]A6583, and (b) [S o]AA6717,6731/[N 
n]A6583 versus [O l]A6300/[N ll]A6583. The asterisks, pluses, and "Y" signs denote the Sis, the SI. 5s, and the S2s of the sample of Nagao 
et al. (2001c), respectively. The data are not corrected for the dust extinction. The data points will move on the diagrams as shown by the 
arrow if the extinction correction of Ay = 1.0 mag is applied. Note that the effect of the extinction correction in the diagram (b) is negligibly 
small. 



As mentioned in the last section, it seems desirable to 
use only forbidden lines for studies of the gas metallicity 
of NLRs in Sis. Furthermore, it seems ideal to use only 
low-ionization forbidden emission lines to discuss the gas 
properties. This is because low-ionization line-emitting re- 
gions and high-ionization line-emitting regions may be spa- 
tially segregated (e.g., Baker 1997; Murayama & Taniguchi 
1998a, 1998b; Nagao, Taniguchi, & Murayama 2000; Na- 
gao, Murayama, & Taniguchi 2001b, 2001c; see also Hes, 
Barthcl, & Fosbury 1993). Taking the above two con- 
straints into account, we examine whether or not the for- 
bidden emission- line flux ratios of [O i]A6300/[N n]A6583, 
[O n]A3727/[N n]A6583, and [S n]AA6717,6731/[N n]A6583 
can be useful estimators for the nitrogen abundance of gas 
clouds in the NLRs. 

In Figure 5, we show the relationship among the three 
forbidden emission-line flux ratios. Here we use the data 
of Nagao et al. (2001c) in order to find possible correla- 
tions among the three flux ratios. Since the data are not 
corrected for the dust reddening, we also show the effect of 
the reddening correction, which is calculated by adopting 
the extinction curve of Cardelli et al. (1989), in Figure 5. 
As shown apparently, positive correlations are seen among 
the three flux ratios. The corresponding correlation co- 
efficients are 0.674 and 0.561 for [O n]A3727/[N n]A6583 
versus [O i]A6300/[N n]A6583 and [S n]AA6717,6731/[N 
n]A6583 versus [O i]A6300/[N n]A6583, respectively. These 
correlation coefficients imply that there are meaningful 
correlations among the three flux ratios. What makes 
these correlations? Or, can these correlations be inter- 
preted as sequences of the nitrogen abundance? In order to 
investigate these issues, we perform photoionization model 
calculations. The model methods and the results are pre- 



sented in the following sections. 

4.2. Photoionization Models 

To investigate the origin of the correlations presented in 
the last section, we carry out photoionization model cal- 
culations by using the publicly available code Cloudy ver- 
sion 94.00 (Ferland 1997, 2000). Here we assume uniform 
density gas clouds with a plane-parallel geometry. The pa- 
rameters for the calculations are (I) the hydrogen density 
of a cloud (nn), (II) the ionization parameter (U), (III) the 
chemical composition of the gas, and (IV) the shape of the 
spectral energy distribution (SED) of the input continuum 
radiation. 

We perform several model runs covering the following 
ranges of parameters: 10 2 cm~ 3 < njj < 10 5 cm~ 3 and 

1Q -3.5 < u < 1Q-2.0 

For the chemical composition of the gas clouds, we as- 
sume the case that the metals are all scaled keeping solar 
proportions except for nitrogen; the nitrogen abundance 
scales with Z 2 (this assumption is altered in some calcula- 
tions presented in §5). We calculate the models covering a 
metallicity range of 0.25 < Z/Z Q < 3.0, which corresponds 
to 8.27 < 12 + log (O/H) < 9.35 and 6.77 < 12 + log 
(N/H) < 8.92. The adopted elemental abundances of the 
solar ones relative to hydrogen are taken from Grevesse 
& Anders (1989) with extensions by Grevesse & Noels 
(1993). We also examine the case that a part of heavy 
elements is depleted into dust grains. In this case, 90% 
of Mg, Si, and Fe, 50% of C and O, and 25% of N and 
S are locked into dust grains, as estimated for the Orion 
H II region (e.g., Baldwin et al. 1991, 1996). This cor- 
responds to the case that the gas clouds in NLRs contain 
unnegligible dust grains. However, we ignore grain opacity, 
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Fig. 6. — Diagnostic diagrams in which the data of Nagao et al. (2001c) are compared with the results of photoionization models. The 
symbols (green) and the arrow are the same as those in Figure 5. The red and the blue lines denote the predictions of the photoionization 
models with the NLSl-like SED and the BLSl-likc SED, respectively. The results of the photoionization models with the non-dcplctcd 
abundances and the depleted abundances are plotted in the left and the right panels, respectively. 



heating and cooling, since their properties in photoionizcd 
gases are highly uncertain. Although the heating process 
by grain photoelectrons is important for enhancing some 
high-ionization emission lines (e.g., Shields & Kennicutt 
1995), it does not affect our results significantly since we 
use only low-ionization emission lines. 

We adopt two types of SED for the calculations; for both 
the SEDs, we adopt the following function: 

'--"--(-i^M-^)^ « 3 > 

(see Ferland 1997). We adopt the following parameter set 
(see Nagao et al. 2001a for more details): (i) the infrared 
cutoff of the big blue bump component, fcTiR = 0.01 Ryd, 
(ii) the slope of the low-energy side of the big blue bump, 
a uv = -0.5, (iii) the UV-to-X-ray spectral index, a ox = 
-1.35, (iv) the slope of the X-ray power-law continuum, a x 
= -0.85 or -1.15, and (v) the characteristic temperature 
of the big blue bump, T BB = 4.9 x 10 5 K or 11.8 x 10 5 
K. Hereafter the continua with (a x , Tbb) = (-0.85, 4.9 x 
10 5 K) and (-1.15, 11.8 x 10 5 K) are denoted as "BLS1- 
like SED" and "NLSl-like SED", respectively. Note that 
the parameter a in equation (3) is determined from the 
adopted value of a ox . The last term in equation (3) is not 
extrapolated below 1.36 eV or above 100 keV. Below 1.36 
eV the last term is simply set to zero. Above 100 keV the 
continuum is assumed to fall off as v~ 3 . 

The calculations are stopped when the temperature falls 
to 3000 K, below which the gas does not contribute signif- 
icantly to the observed optical emission-line spectra. 

4.3. Results of Model Calculations 



In Figure 6, we present the results of the models with 
the mctallicity of Z = Zq. These model results are used 
to investigate whether or not the gas density, the ioniza- 
tion parameter, or the combination of both two parameters 
drives the correlations presented in §4.1. Here the models 
with the NLSl-like SED and the BLSl-like SED are pre- 
sented. It is apparently shown that the emission-line flux 
ratios of [O i]A6300/[N n]A6583, [O n]A3727/[N n]A6583, 
and [S n]AA6717,6731/[N n]A6583 are insensitive to the 
ionization parameter for any type of the adopted ioniz- 
ing SEDs. The predicted emission- line flux ratios vary by 
less than a factor of 2 when the ionization parameter is 
changed between 10~ 35 and 10~ 20 . This is mainly be- 
cause we use only low-ionization forbidden emission lines. 
Although the three emission-line flux ratios are rather sen- 
sitive to the gas density, the model grids for varying gas 
densities are perpendicular to the observed correlations for 
any type of the adopted ionizing SEDs, as shown in Fig- 
ure 6. These results do not depend on whether or not the 
heavy elements are depleted into dust grains, as shown in 
Figure 6. Therefore, we conclude that the observed cor- 
relations among the three emission-line flux ratios are not 
significantly attributed to the effects of the gas density, 
the ionization parameter, or both. 

Then, are the correlations attributed to the nitrogen 
abundance (i.e., the gas metallicity) ? In Figure 7, we 
show the results of model calculations in which the gas 
metallicity is varied between 0.25 < Z/Zq < 3.0. Here 
the NLSl-like SED and the BLSl-like SED are adopted 
for the ionizing continuum together with U = 10~ 30 . It 
is shown that the metallicity sequences can well explain 
the observed correlations in the flux ratios of the low- 
ionization forbidden lines for both cases, the models with 
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Fig. 7. — Same as Figure 6, but other model grids are overdrawn. As for the photoionization models, the ionization parameter is fixed to 
be U = 10 — 3 . The metallicity is varied between 0.25 < Z/Zq < 3.0. The symbols, the lines, and the arrow are the same as those in Figure 
6. 



the non-depleted abundances and with the depleted abun- 
dances. We thus conclude that the observed correlations 
among the flux ratios of the low-ionization forbidden lines 
are mainly attributed to the variety of the gas metallic- 
ity. Since the narrow component of permitted lines is not 
used and only forbidden lines are used, these diagrams pre- 
sented in Figure 7 can be powerful diagnostics for the gas 
metallicity in the NLR of AGNs, especially for Sis. 

Here we mention that the flux ratio of [O i]A6300/[N 
n]A6583 seems to be the most appropriate diagnostic for 
gas metallicity among the three emission-line flux ratios. 
This is because the flux ratio of [O ll]A3727/[N n]A6583 
and [S il]AA6717,6731/[N li]A6583 are sensitive to the gas 
density while that of [O i]A6300/[N ll]A6583 is not, unless 
the gas density is high (> 10 4 cm -3 ). The flux ratio of 
[O ll]A3727/[N n]A6583 is sensitive to the dust extinction, 
additionally. 

4.4. Application to NLRs in NLSls and BLSls 

Now we apply our new method to the investigation 
of the gas metallicity of NLRs in NLSls and BLSls. 
In Figure 8, we plot the data of Rodriguez- Ardila et 
al. (2000) on the diagnostic diagrams presented in §4.3. 
It appears that the NLSls have smaller flux ratios of 
[O i]A6300/[N n]A6583, [O ll]A3727/[N n]A6583, and [S 
n]AA6717,6731/[N n]A6583 than the BLSls. The average 
and median values of these three emission-line flux ratios 
for the NLSls and BLSls are given in Table 1. These 
values imply that there are differences in the emission- 
line flux ratios between the NLSls and the BLSls, though 
the statistical significance is low; i.e., Pks = 0.133, 0.533, 
and 0.767 for the flux ratios of [O i]A6300/[N n]A6583, [O 
n]A3727/[N n]A6583, and [S n]AA6717,6731/[N n]A6583, 
respectively (see Table 1). Here it should be noted that 



the photoionization models predict larger ratios of [O 
l]A6300/[N n]A6583 for NLSls than for BLSls, if the dif- 
ference in the SED of ionizing photons between NLSls and 
BLSls is taken into account. The data of Rodriguez- Ardila 
et al. (2000) exhibit the opposite trend from the predic- 
tion of the photoionization models. It is thus suggested 
that the gas metallicity in NLRs of the NLSls is possibly 
higher than that of the BLSls, at least for the sample of 
Rodriguez- Ardila et al. (2000). 

In Figure 9, we show the relationship between the 
FWHM of the broad component of the Ha emission and 
the flux ratio of [O i]A6300/[N n]A6583 for the sample 
of Rodriguez- Ardila et al. (2000). Being different from 
Figures 2 and 3, it appears to be a positive correlation be- 
tween them. Its correlation coefficient is 0.681, which im- 
plies that there is a meaningful correlation between them. 
This is consistent with the result of Wills et al. (1999); 
that is, the object with narrower BLR emission tends to 
exhibit higher nitrogen abundance. 

5. DISCUSSION 

As investigated in §3 and §4, the inferred difference in 
the gas metallicity of NLRs between NLSls and BLSls 
depends on the adopted diagnostics. When equation (1) 
is used to estimate the gas metallicity, it is suggested 
that NLSls tend to have metal-poorer NLR gas than 
BLSls, though the difference is not significant. Con- 
trary to this result, the diagnostics examined in §4 suggest 
that NLSls appear to possess metal- richer NLR gas than 
BLSls. There appears to be no significant difference in 
the estimated metallicity by using equation (2). Why the 
results are inconsistent in each other? Which is a reliable 
diagnostics to estimate the gas metallicity in NLRs? Here 
we discuss these issues in order to understand how the gas 
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Fig. 8. — Same as Figure 7, but the data of Rodriguez-Ardila et al. (2000) are plotted instead of the data of Nagao et al. (2001c). The 
lines and the arrow are the same as those in Figure 6. The red circles and the blue triangles denote the data of the NLSls and the BLSls in 
the sample of Rodriguez-Ardila et al. (2000), respectively. 



1 



So. 



1 cr 10 
FWHM(Ha broad ) 



Fig. 9. — Diagram of the flux ratio of [O l]A6300/[N n]A6583 versus the FWHM of the broad component of the Ho emission. The symbols 
are the same as those in Figure 2. 
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metallicity is different, or not, between NLSls and BLSls. 

As discussed by Nagao et al. (2001c), it is sometimes 
difficult to measure fluxes of the narrow component of 
permitted lines for Sis. Since the method using cither 
equation (1) or (2) requires the flux of the narrow compo- 
nent of permitted lines, the oxygen abundances derived at 
§3 can suffer from this effect and thus he derived values 
may contain large errors. In addition to this uncertainty, 
there is another effect which may diminish the reliabil- 
ity of the metallicity diagnostics adopted in §3. Veron- 
Cetty et al. (2001) reported that the broad Balmcr lines 
in the spectra of NLSls are well fitted by a Lorentzian 
profile rather than a Gaussian profile (see also Moran et 
al. 1996; Leighly 1999b; Goncalves et al. 1999; Sulentic 
et al. 2002). When the broad component is fitted by a 
Lorentzian profile instead of a Gaussian profile, the mea- 
surement of the narrow component of the line is affected. 
This results in larger ratios of [N ii]A6583/Ha narrow and [O 
in]A5007/H/3 narrow compared to the case when the broad 
component of permitted lines is fitted by a Gaussian pro- 
file as done by Rodriguez-Ardila et al. (2000) (see Figure 
8 of Veron-Cetty et al. 2001). Therefore, the ratios of [N 
ii]A6583/Ha narr0 w and [O iii]A5007/H/3 narrow in the data 
of Rodriguez- Ardila et al. (2000) may be underestimated 
if a Lorentzian profile is a better function to measure the 
broad component of permitted lines than a Gaussian pro- 
file which Rodriguez-Ardila ct al. (2000) adopted. In 
order to see how the underestimation of the flux ratios 
of [N n]A6583/Ha narrow and [O m]A5007/H/3 narrow affects 
the inferred oxygen abundance, we consider the case in 
which the two flux ratios are underestimated by a fac- 
tor of 10. We assume [N n]A6583/Ha narrow = 2.0, [O 
m]A5007/H/3 narrow = 10.0, and [O n]A3727/[0 m]A5007 
= 0.1, which would be measured as [N n]A6583/Ha narrow 
= 0.2, [O m]A5007/H/3 narrow = 1.0, and [O n]A3727/[0 
in]A5007 = 0.1, respectively. In this situation, the oxy- 
gen abundances are also underestimated by a factor of 
-0.28 for log (0/H)sscki and by a factor of -0.30 for 
log (0/H)sscK2- This result suggests that we might un- 
derestimate the oxygen abundance in §3 if the Lorentzian 
profile is the representative function for the broad compo- 
nent of permitted lines. This may be the reason why the 
results derived in §3 and §4 are inconsistent; the oxygen 
abundance of gas in NLSls may be selectively underes- 
timated when the SSCK method is applied, due to the 
Lorentzian profile of broad components of permitted lines 
in the spectra of NLSls. Although the reason why the 
broad component of the permitted lines is well fitted by 
Lorentzian is not clear, Dumont & Collin-Souffrin (1990) 
showed that the Lorentzian profile could be produced if the 
BLR emission is produced at gas clouds in a disk-like ge- 
ometry around a supermassive black hole (see also Veron- 
Cetty et al. 2001; Sulentic et al. 2002). We do not discuss 
the origin of the Lorentzian profile further since this issue 
is beyond the purpose of this paper. Since the diagnostics 
investigated in §4 consist of only forbidden lines and thus 
are free from the above uncertainties, we conclude that the 
results obtained in §4 are more reliable than those obtained 
in §3. Therefore, it can be concluded that the metallicity 
of gas clouds in NLRs is higher in NLSls than in BLSls. 

Here we mention that the diagnostics investigated in §4 
may also have a possible drawback. As shown in Figure 
7, the correlations among the emission-line flux ratios of 



[O i]A6300/[N n]A6583, [O n]A3727/[N n]A6583, and [S 
n]AA6717,6731/[N n]A6583 can be explained by the pho- 
toionization models in the parameter ranges of «h ~ 10 2 ~ 5 
cm^ 3 and 0.5 < Z/Zq <J 3.0. The inferred range of the gas 
density is consistent with the previously reported values 
(e.g., Tsvetanov & Walsh 1992; Storchi-Bergmann, Wil- 
son, & Baldwin 1992; Schmitt et al. 1994). However, the 
model values of metallicity would contain possible uncer- 
tainties if the nitrogen production is delayed relative to the 
oxygen production. The delay is significant if the timescale 
for the increase in the metallicity by star formation activ- 
ity is comparable to or shorter than the lifetime of the 
nitrogen-producing stars, i.e., ~ 10 s yr (e.g., Hamann & 
Ferland 1993, 1999; Henry, Edmunds, & Koppen 2000; see 
also Hamann ct al. 2002). This effect results in that the 
abundance ratio of N/H become to be below the solar value 
at Z = Z Q . Although this offset is only weakly present in 
the data of H n region (van Zee et al. 1998), it could be 
substantial in AGNs (see, e.g., Hamann & Ferland 1999). 
Taking this effect into account, one can formulate the ni- 
trogen abundance of the gas as 

log jmk = 2log{z/z&) - q (4) 

where q is the logarithmic offset. Hamann et al. (2002) 
mentioned that the q value might reach up to —0.2-0.5 
for gas in AGNs. In order to investigate the effect of this 
delay on our diagnostics, we perform further model calcu- 
lations. The adopted parameters are the same as those of 
the models presented in Figures 7 and 8 except for the q 
value, which is altered from q = to q = 0.5. The results 
are shown in Figure 10. The diagrams in Figure 10 suggest 
that the gas metallicity of NLRs is 1.0 <; Z/Z Q <; 3.0 for 
the BLSls while Z/Z Q £ 2.5 for the NLSls. These values 
are larger than those suggested by the diagrams in Figure 
8. Therefore, the diagnostics investigated in §4 may be 
invalid when we estimate the absolute values of the gas 
metallicity since the true q value for the NLR in AGNs 
is unclear. We can, however, use the diagnostics to com- 
pare the gas mctallicities between NLSls and BLSls unless 
the q value is systematically different between NLSls and 
BLSls. Since there is no reason to exist such a differ- 
ence in the q value between the two populations of Sis, 
the conclusion that the gas metallicity in NLRs tends to 
be relatively higher in NLSls than in BLSls seems valid, 
even if the uncertainty of the q value is taken into account. 

Finally, we discuss the effect of the inhomogeneity of 
NLRs on our results. It is known that there are significant 
density gradients in NLRs (e.g., Kracmcr ct al. 2000). 
Since the [O i]A6300 transition has higher critical den- 
sity (1.8 x 10 6 cm" 3 ) than the transitions of [O n]A3727, 
[N n]A6583, [S n]A6717, and [S n]A6731 (4.5 x 10 3 cm" 3 , 
8.7 x 10 4 cm~ 3 , 1.5 x 10 3 cm~ 3 , and 3.9 x 10 3 cm~ 3 , respec- 
tively), gas clouds with ran > 10 5 cm~ 3 may selectively 
radiate the [O i]A6300 emission. Therefore our one-zone 
photoionization models may be inappropriate to investi- 
gate the gas metallicity of NLRs if such a high-density gas 
cloud contributes to emission-line spectra of NLSls and 
BLSls differently. Indeed there are significant differences 
in the emission-line fluxes of some high-ionization emis- 
sion line such as [Fe vn]A6087 and [Fe x]A6374 between 
Sis and S2s, which is due to the different contribution of 
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Fig. 10. — Same as Figure 9, but the models with q = 0.5 are drawn instead of the models with q = 0. 



dense gas clouds in NLRs (e.g., Murayama & Taniguchi 
1998a; Nagao et al. 2000, 2001c). However, it has been al- 
ready shown by Nagao et al. (2000, 2001c) that there is no 
systematic difference in the high-ionization emission lines 
between NLSls and BLSls. This suggests that the contri- 
bution from the dense gas clouds to emission-line spectra 
is not so different between NLSls and BLSls. Therefore, 
our conclusions seem valid even if possible density inho- 
mogeneity is taken into account. 

6. SUMMARY 

In order to explore possible differences in the gas metal- 
licity between NLSls and BLSls, we have focused on the 
gas in NLRs. It is found that the results depend on the 
applied diagnostics. The emission-line flux ratio of [N 
ii]A6538/Ha narrow , which has been often used to discuss 
the gas metallicity of NLRs in S2s, is not significantly dif- 
ferent between NLSls and BLSls. Some other diagnostics 
in which the narrow component of permitted lines is used 
also suggest that there is little or no systematic difference 
in the gas metallicity between NLSls and BLSls. On the 
other hand, the diagnostics which consist of only forbid- 



den emission lines, which are newly proposed in this paper, 
suggest that the gas metallicity of NLRs in NLSls tend to 
be higher than that in BLSls. This complex situation may 
be solved if the Lorentzian profile, rather than the Gaus- 
sian profile, is the representative function for the broad 
component of permitted lines in the spectra of NLSls. 
Since the results derived by using only forbidden emission- 
line fluxes do not suffer from the uncertainty of the fitting 
function for the broad component of permitted lines, it 
can be concluded that the gas metallicity of NLRs tends 
to be higher in NLSls than in BLSls. 
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